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Abstract 

An approach to calculating approximate solutions to the continuum Schwinger- 
Dyson equations is outlined, with examples for 4 in D = 1. This approach is 
based on the source Galerkin methods developed by Garcia, Guralnik and Law- 
son. Numerical issues and opportunities for future calculations are also discussed 
briefly. 

1 Introduction 

The now-conventional technique of numerical calculation of quantum field theory in- 
volves evaluating the path integral 



& 1 ! Z{J) = I d<f> exp [-5(0) + J<f>] (1) 



X 



Jr 

on a spacetime lattice using Monte Carlo integration methods. Monte Carlo methods 
have been successful for an interesting class of problems; however, these techniques 
have had less success in evaluating theories with actions that are not manifestly posi- 
tive definite or which have important effects from the details of fermionic interactions 
beyond the quenched approximation. In these cases, we have either a basic algorithmic 
difficulty or an inadequacy of compute power. 

We are developing an alternative computational method which works both on the 
lattice and the continuum and which handles fermions as easily as bosons. Furthermore, 
our "source Galerkin" method is less restrictive as to the class of allowed actions. Source 
Galerkin tends to use significantly less compute time than Monte Carlo methods but can 
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consume significantly larger amounts of memory.^ This talk is confined to continuum 
applications; examples of lattice calculations have been given elsewhere |7|, . While 
it is not clear that Source Galerkin can replace Monte Carlo techniques, it appears that 
it will be able to solve some problems which are currently inaccessible. 

Our approach begins with the differential equations satisfied by the vacuum func- 
tional Z for a quantum field theory with external sources. For the sake of simplicity, in 
this talk we will mostly confine our attention to (f) 4 interactions. To date, we have also 
studied non-linear sigma models and four-fermion interactions and have gauge theory 
calculations in progress. 

The vacuum persistence function Z for a scalar field <fi with interaction g0 4 /4 cou- 
pled to a scalar source J{x) satisfies the equation: 

The source Galerkin technique is designed to directly solve functional differential equa- 
tions of this type. Before we proceed to outline a solution technique, it is essential to 
point out that this equation by itself does not uniquely specify a theory ||. This is 
dramatically illustrated by considering the special case of the above equation limited 
to one degree of freedom (zero dimensions). 

, r9 dZ d 3 Z T 

M u + ^- JZ = < 3 > 

This is a third order differential equation, and therefore possesses three independent 
solutions. It is easy to see that one of the solutions for small g asymptotically approaches 
the perturbation theory solution, while the second asymptotically approaches the loop 
expansion "symmetry breaking" solution for small g and the third solution has an 



This can happen because the iterative process for improving a source Galerkin calculation involves 
successively higher-point Green functions, whereas calculation of these correlations in a Monte Carlo 
scenario is optional. 
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essential singularity as g becomes small. The situation becomes much more interesting 
for finite dimensions where the infinite class of solutions coalesces or becomes irrelevant 
in a way which builds the phase structure of the field theory. Any numerical study must 
be cognizant of the particular boundary conditions and hence solution of the class of 
solutions that is being studied. Care must be taken to stay on the same solution as 
any iterative technique is applied. The solutions discussed in this talk will be the usual 
ones which correspond to the symmetry preserving solutions obtained from evaluating 
a path integral with real axis definitions for the regions of integration. These are the 
solutions that are regular in the coupling, g, as it approaches zero. 
We, for convenience, write (|2]) in the form: 



The source Galerkin method is defined by picking an approximation Z*(j) to the solu- 
tion Z(j) such that 



where R is a residual dependent on j. We pick the parameters of our approximation to 
make this residual as small as possible on the average. To give this statement a mean- 
ing, we define an inner product over the domain of j: i.e. (A,B) = j dfi(j) A(j)B(j). 
In addition we assume we have a collection of test functions which are members of a 
complete set) {(fi{j)} The source Galerkin minimization of the residual R is imple- 
mented by setting the parameters of our test function Z*(j) so that projections of test 
functions against the residual vanish so that | \Z* — Z\ [2 — > as the number of test 
functions — > 00. 

The equations defining the quantum field theory are differential equations in the 
field sources and spacetime. While it is straightforward to deal with the spacetime 
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problem by resorting to a lattice, we can remain in a continuum formulation by tak- 
ing advantage if our knowledge of functional integration. We know how to evaluate 
Gaussian functional integrals on the continuum: 



exp 



j(x)A(x,y)j(y)+ / j(x)P(x 



exp 



VdetA 

Consequently, we can evaluate integrals of the form: 

[dj]exp [-f(x)/e 2 ] P{ 3 ) 



(3*(x)A- 1 (x,y)(3(y) 



(6) 



(7) 



Using this we can define an inner product of sources on the continuum as follows: 



(j(xi) ■ ■ ■ j{x n ), j(y 1 ) ■ ■ -j{y m )) 3 



e n+m 5 + {xi ■ ■ ■ x n yi ■ ■ - y m } n + m even 







otherwise 



where we have absorbed a factor of 2 by redefining e. 5 + is defined by 

5 + {xa(3 • • • } = S(x — a)5 + {/3 • • • } + S(x — [3)5 + {a •••} + ••• 

5 + {xa} = 6(x — a). 



(9) 
(10) 



In addition to this inner product definition, we need good guesses for approximate 
form for Z* and numerical tools to calculate, symbolically or numerically, various func- 
tions and their integrals, derivatives, and so on. For most of our calculations we have 
found it very useful to choose a lesser known class of functions, with very suitable 
properties for numerical calculation, known as Sine functions. We take our notation 
for the Sine functions from Stenger ||: 

sin(7r(x — kh)/h) 



S(k,h)(x) 



tt(x — kh)/h 



(11) 
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Some of the identities that Sine approximations satisfy are given below: 

S(k, h)(lh) = S k i (12) 
/ S(k,h)(x)S(l,h)(x)=5 kl (13) 

J X 

r - 

/ F{x) ~hJ2 F ( kh ) ( 14 ) 

^ x k=-N 
N 

F{x)tt F{kh)S{k,h){x) (15) 

k=-N 
N 

F'(x)& F{kh)S'{k,h){x) (16) 

k=-N 

N 

Fl»\lh) « F (kh)6™ (17) 

k=-N 

These properties, which are proven and expanded upon greatly in ||, make these 
functions very easy to use for Galerkin methods, collocation, integration by parts, and 
integral equations. 

With the definition of a norm and set of expansion functions, we can postulate an 
ansatz for Z 



Z* = exp 



^2 / j( x ) G 2( x ~y)j(y) + 

J xv 



' xy 

where each Green function, G n , is represented by a ci-dimensional Sine expansion. For 
d = 4: 

G 2 (x-y) = J2Gf l 

S(i, h)(x° - y°)S(j, h){x l - y x )S{k, h)(x 2 - y 2 )S(l, h)(x 3 - y 3 ) (19) 

It is easy to examine this expansion in the case of a free field where we limit the 
approximation to the terms quadratic in j(x). While the example is trivial, it shows 




Figure 1: Mass, m*, versus distance, x, for two-dimensional free scalar field. Note 
breakdown near origin (approximation of ^-function) and at large distance (spatial 
truncation). m = 0.5. 

that, as always, numerical approximations must be handled with care. Results of this 
calculation are shown in Figure [I]. 

This straightforward expansion works for interacting theories but with more than 
G2, the computational costs and storage costs become overwhelming: G2n requires 
jy(2n-i)d s t ora g e un fts. There are many ways that storage costs can be reduced, but in 
general these approaches are difficult, not particularly elegant, and eventually reaches a 
limit due to the exponential growth in the number of coefficients of the representation. 

We can use our knowledge of the spectral representations of field theory and graphi- 
cal approaches to introduce a much more beautiful and intuitive approach to producing 
candidates for Z* . We introduce regulated Lehmann representations. These build in 
the appropriate spacetime Lorentz structure into our approximations and remove the 
growth of operational cost with spacetime dimension shown by our previous naive de- 
composition into complete sets of functions. Any exact two-point function can be 
represented as a sum over free two-point functions. We choose as the basis of our 
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Exact 


0.654038297612956387655447 


N = 10 
20 
40 
80 
100 


0.686283404027993439164852 
0.Q55465233659481 763033090 
0.65403 7798068595478871392 
0.654038297613089305988864 
0.654038297612956949544002 



Table 1: Convergence of Sine approximation to integral. (Exact integral calculated 
using Maple V, with 40 digits of precision.) x 2 = 10, m = 1, A 2 = 10. 



numerical solutions, a regulated Euclidean propagator structure: 

ip-x— p 2 /A 2 

'(dp) I dse ip - x - pI/A2 - s{p2+m2) 
Jo 

d/2 



A(m; x) = I (dp) 



p2 _|_ m 2 



ds 



71 



s + 1/A 2 



"4(s + l/A^) 



(20) 



(21) 



(2tt) - jo 

This regulation assures that we never have to deal with infinities in any calculated 
amplitude as long as we keep the cutoff finite. 

This integral can be approximated using Sine methods 



A(m; x) 



h 



N 1 

-V — 

d / j p kh 



(2,71 J" " e 
v I k=-N 



71 


d/2 


exp 


[z k + l/A\ 



X 



-z k m - 



4(z k + 1/A 2 ) 



:22) 



The example in Table |IJ demonstrates that we can have as many digits as are necessary 
for the calculation, with the associated increase in compute time. For practical pur- 
poses, 80 terms is appropriate for most hardware floating-point representations. Thus 
we have a form for a two-point scalar Green function, regulated by the scale A 2 with 
constant computational cost regardless of spacetime dimension. We can take derivatives 
explicitly or by construction: 



d 2 A(m; x) = m 2 A(m; x) — S(x) 



(23) 



where S(x) = e 



-x 2 A 2 /4 
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Figure 2: One dimensional 4 mass gap versus coupling (dashed line gives exact from 
Hioe and Montroll, 1975) 



From this representation, we can directly construct a fermion two-point function: 



These representations mean that free scalar and free fermion results are exact and 
immediate in any Galerkin evaluation of these trivial cases. Furthermore, because of 
this simplicity, we have the basis for a complete numerical approach to conventional 
perturbation theory. 

2 Results with Lehmann representation: </> 4 

We itemize some results obtained using a regulated single propagator with parameters 
set by the Source Galerkin method. At lowest order, our ansatz for the generating 
functional is 



Results for this ansatz are given in Figure ||. These results are strikingly accurate and 




(24) 




(25) 
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Figure 3: Additional connector-based ansatzen for the four-point function, H, in A0 4 . 
In the bottom row, we have two contact ansatzen on the left, followed by two mediated 
ansatzen. 

can matched up essentially exactly with results of Monte Carlo calculations in two and 
higher dimensions. 

We can enhance these results by including additional 4 source terms in Z*. Some 
simple additional terms that we include with weights and masses to be calculated using 
the Source Galerkin technique are the terms of the forms given in Figure |3|. The effect 
of adding a fourth order term is shown in Figure f|. 

In addition to the illustrations given here, we have examined {ipip) 2 in the mean 
field in two dimensions and have found rapid convergence to the known results from 
large- N expansions. It therefore appears that, at least for the simple cases studied, we 
have produced a numerical method which draws on the structural information already 
known in general through symmetry and spectral representations which when combined 
with Galerkin averaging to set parameters converges with very simple guesses for the 
vacuum amplitude to known correct answers produced through other methods of so- 
lution including Monte Carlo methods. More complicated gauge problems are under 
study. 
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2 4 

Coupling, A 



Figure 4: Comparison of fom-H approximation with lowest order and exact answer, in 
one-dimension (A 2 = 70). 
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3 Numerical issues 



In this very brief presentation we have avoided discussion of many of the difficult 
numerical issues involved in constructing this approach. We note some of these issues 
here without discussion to have on record: 

• Both interpolative and spectral problems result in medium- to large-scale non- 
linear systems; systems solvable using many variable Newton's method 

• Finite storage is the key constraint for interpolative representations, which must 
be constrained to two-point "connectors" , particularly in high dimensions 

• Storage is a non-issue for spectral representations (memory use entirely for caching 
calculated quantities) 

• Both methods are also time-bound to connector-based representations for higher 
point functions 

• Time cost from internal loops; however, algorithm can be made parallel via par- 
titioning of sums 

• Resolution of elementary pole structures (i.e. differentiating between S and (p 2 + 
m 2 ) -1 ) may be addressable with arbitrary precision libraries 

• Arbitrary precision may also be useful for regulated perturbation calculations. 

These numerical issues, and related general numerical techniques for source Galerkin 
are discussed in greater detail in fJXJ, |l| f|. 
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4 Conclusions 



We have discussed a method for numerical calculations for field theories on the contin- 
uum; this method being based on the source Galerkin technique introduced in . The 
direct approach, using Sine functions for interpolation, is effective but will ultimately 
be limited by the finite nature of current computational resources. The Lorentz in- 
variant regulated representation derived from Lehmann representation does not suffer 
from these limitations, and is applicable to both perturbation and mean-field-based 
ansatzen. This approach has the computational advantages of minimal memory uti- 
lization and parallelizable algorithms and also allows direct representation of fermionic 
Green functions. Finally, a number of useful peripheral calculations can made using 
this appoximate representation: one can calculate diagrams in a regulated perturbation 
theory, as well as calculating dimensionally regularized loops numerically. In general, 
this technique of evaluating field theories takes advantage of the symmetries of the 
Lorentz group; future work includes the extension of the method to more general in- 
ternal groups, such as gauge groups or supersymmetry. 
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